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This thesis designs and implements the Coordipated falaag SE SMATESCHOOL 
Rescue (SAR) System (COINSS). This system provides several important features not 
provided by current inland SAR computer systems. First is the ability to model 
movement of the target. Second is modeling the effect terrain has on the movement of 
the target. Third is the visual presentation of a probability map, a color display showing 
the probability that the target is located at various geographic positions. COINSS is 
developed in the Java programming language. It is designed to be implemented with a 
map-based planning system using loosely coupled components. COINSS provides the 
initialization, movement, and search algorithms which are used by the planning system to 
support the search operation. The initialization algorithms define the search area where 
the SAR operation will occur. Initial areas are defined for the target. COINSS models 
the movement of the target as a discrete trme Markov chain. Bayes theorem 1s used to 
update the probability map when negative search information is provided. This thesis 


will improve inland SAR operations by providing the first model with an interactive 


graphical user interface and a model of target movement. 
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EXECUTIVE SUMMARY 


The National Search and Rescue Plan establishes the United States Air Force as 
the executive agent for Inland Area Search and Rescue (SAR). This area encompasses 
the Continental United States, except Alaska, and waters under the jurisdiction of the 
United States. The Air Force established the Air Force Rescue Coordination Center 


(AFRCC) as the contact point for federal resources for inland SAR [18]. 


In addition, the National Search and Rescue Plan dictates that Search and Rescue 
Coordinators should provide for the fullest cooperation and use of additional Federal 
military, civil, State, local, and private agencies as may be necessary and practicable [18]. 
The On-Scene Commander is responsible for overseeing the execution of the entire 
search plan. The Commander uses search tables to formulate each assignment in the 
overall search plan. The individual search assignments are then provided to each search 
unit for implementation. The results of the individual search assignments are 
communicated back to the Commander and used to update the search scenario. Based 
upon the updated scenario, the Commander reevaluates the search plan, creates new 


search assignments, and the process is repeated. 


Currently, numerous diverse software tools and spreadsheets are available to help 
do numeric planning calculations and table referencing. These stand-alone aids compute 
the data for the Commander, who must then interpret the data and incorporate it into the 
search plan. Additionally, management programs are available for ongoing analysis of the 


search plan. 


1X 


However, despite the availability of a variety of search software programs, three 
key factors prevent their widespread use and limit their usefulness during an active inland 


search operation. These factors are: 
1. the lack of a graphical mapping interface, 


2. the lack of inter-operability among the computer systems of cooperating 


agencies, and 
3. the inability to account for a mobile target. 


This thesis designs and implements the Coordinated Inland Area Search and 
Rescue (SAR) System (COINSS). This system accounts for the three factors not 
provided by current inland SAR computer systems. First, COINSS provides a color 
display of the target probability distribution, which can be used as an overlay on a map- 
based planning system. Second, it is developed in the Java programming language and is 
designed to be implemented with a map-based planning system using loosely coupled 
components. Finally, it models the movement of the target as a discrete time Markov 
chain. 

COINSS provides the initialization and search algorithms which are used by the 
planning system to support the search operation. The initialization algorithms define the 
search area where the SAR operation will occur. Initial areas are defined for the target. 
Bayes theorem is used to update the probability map when negative search information is 
provided. 

This model will improve inland SAR operations by providing the first model with 


an interactive graphical user interface and a model of target movement. 
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THESIS DISCLAIMER 


The reader is cautioned that the computer programs developed in this research 
may not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and logic 
errors, they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 


X111 








I. INTRODUCTION 
A. BACKGROUND 
1. The Responsibility of Land-Based Search and Rescue (SAR) 


The National Search and Rescue Plan establishes the United States Air Force as 
the executive agent for Inland Area Search and Rescue (SAR). This area encompasses 
the Continental United States, except Alaska, and waters under the jurisdiction of the 
United States. The Air Force established the Air Force Rescue Coordination Center 


(AFRCC) as the contact point for federal resources for inland SAR [18]. 


In addition, the National Search and Rescue Plan dictates that Search and Rescue 
(SAR) Coordinators should provide for the fullest practicable cooperation and use of 
additional Federal military, civil, State, local, and private agencies as may be necessary 
and practicable [18]. To this end, the AFRCC has established agreements with each 
State, which describe the responsibilities for different emergencies. However, these 
agreements are not consistent from State-to-State, and often do not fully address the 


availability of local and private agencies to assist with the search effort. 
2. SAR Tools and Process 


Search planning is a well-researched and documented discipline. The National 
Search and Rescue Plan provides a complete set of tables and charts for creating and 


calculating search areas and search patterns. 


The On Scene Commander is responsible for overseeing the execution of the 


entire search plan. The Commander uses search tables to formulate each different search 


assignment in the overall search plan. The individual search assignments are then 
provided to each search unit for implementation. The results of the individual search 
assignments are communicated back to the Commander and used to update the search 
scenario. Based upon the updated scenario, the Commander reevaluates the search plan, 


creates new Search assignments, and the process is repeated. 


The bulk of search planning is still done on paper charts and maps. The 
Commander must take into account a large range of factors such as the target's last 
known position, time of that position, the terrain that is being searched, the type of target 
that is being sought, the number, experience, and composition of the available search 
units, etc. The Commander uses this information, as well as personal experience and all 
available charts and tables, to lay out the search area, divide the area into appropriate 
search sectors, assign individual search teams to cover the sectors, and establish the 


search pattern(s) that will be used’ by the search teams. 


Currently, numerous diverse software tools and spreadsheets are available to help 
do numeric planning calculations and table referencing. These stand-alone aids compute 
the data for the Commander, who must then interpret the data and incorporate it into the 


search plan. 


Additionally, management programs, such as CASIE III [3], are available for 
ongoing analysis of the search plan. These programs allow the searcher to define a 
search region and divide that region into search areas. Each area is assigned a probability 
that the target is located within that area. The searches are then conducted in the various 


search areas and the results used to update the probabilities by applying Bayes theorem. 


An unsuccessful search in an area reduces the probability that the target is in that area. 
This reduction in probability is directly related to probability of success of the search 


method used. 
B. PROBLEM 


Despite the availability of a variety of search software programs, three key factors 
prevent their widespread use and limit their usefulness during an active inland search 


operation. These factors are: 
1. the lack of a graphical mapping interface, 


2. the lack of inter-operability among the computer systems of cooperating 


agencies, and 
3. the inability to account for a mobile target. 


CASIE III is an example of current software that provides the Commander with 
the raw data required to make search decisions. The program takes user input search 
areas and target location probabilities and the program provides, as output, probabilities 
that the target is located within the predefined search areas. When the results of a search 
are entered, the program then updates the probabilities, and outputs these revised values. 
The Commander is then responsible for translating that information to a physical map 
representation. Then the Commander assesses the search region and evaluates the next 
step to be taken in the search execution. The new search assignments, and eventually 


search results, must then be entered back into the program, where the cycle resumes. 


The second factor that must be considered is the wide assortment of resources that 
can be brought to help execute aSAR case. The range of available personnel resources 
goes from trained, full-time federal and state professionals all the way to sparsely trained 
local volunteer organizations. The computer resources that will be available from these 
sources are unpredictable at best. In order for a system to be useful for all organizations 
participating in the search, the program must be able to run or be viewed on any 


computer system. 


The third factor is the abilities of the program itself. No model currently exists 
which accounts for a moving target on land. Probabilities are assigned to the various 
areas of the search region. After these probabilities are assigned, the target is assumed to 
be immobile. The changes in the probabilities in each area are only the result of 


unsuccessful searches that were conducted. 
C. APPROACH 


This thesis develops the Coordinated Inland Area SAR System (COINSS) to 
provide one method of resolving the problems addressed above. The design of this 


computer program was based upon several key elements. 
1. The Key Elements 
a. Portable 


Due to the urgent and unpredictable nature of SAR, it is not possible to 
prepare in advance for all possible scenarios. The Commander must make best 


use of the resources available at the time. The computer program must be small 


enough to be easily and quickly loaded from a disk or downloaded over a network 
and provide a user interface that requires little user training. Additionally, the 


program should be able to operate on a wide variety of systems. 
b. Graphical interface 


A graphical interface would allow the Commander to do the SAR 
planning on the system, and not require a transfer of information between a paper 
chart and the computer model. The visual display of probabilities would provide 
a better description of the search progress. Additionally, the direct interface 
between the visual display and the search logic would provide a quicker update to 


the search execution once other search information has been received. 
c. Interoperable with available maps 


Again, SAR is unpredictable. We cannot know where the next SAR case 
will occur. The program should allow for use of maps from any electronic source 


to provide the graphical planning interface. 
d. Allow for target movement 


A target may not remain stationary. Ignoring the motion of the target in 
the execution of the SAR plan significantly impacts the location of the search 


areas, the search method used in those areas, and the results of the search. 


e. Common picture view 


One advantage of the portability element is that al] units coordinating on a 
search would be able to view the planning tool to assist in their search execution. 
The Commander would maintain and update the operational view of the system, 


which would be distributed to each search unit. 
2. Two Concepts of Movement 


Two choices are available to describe the target movement through the search 
area. Wagner [19] explains one method as being represented by a bundle of simulated 
tracks that encompass the many possibilities for target initial position and movement. 
The number of necessary tracks depends upon the complexity of the target being 
modeled, and should be set to a sufficiently large number in order to cover the complete 
range of possible initial conditions. A probability 1s associated with each track. This is 
the probability that the target is following the respective track. Wagner further explains 
how the tracks are monitored and updated, based upon negative search results. These 
tracks move over a grid search area, and require the use of random number generation to 


determine the progressing course of the simulated tracks. 


The second method is a pure stochastic representation of movement [19]. A grid 
is placed over the search area. Each grid square, S,,, contains a velocity array, Vx,yij, 
which contain elements that are the probability that the target is located in the square, 
moving with velocity i in the east-west direction, and velocity j in the north-south 


direction. The total probability for each square is computed by: 


lod i evar 
iJ 


During time-based updates, the probabilities are moved from square to square, 
based upon their velocity, and then the velocity state is updated based upon a Markov 
chain. Bayes theorem is used to update the probability distributions when negative 


search information is received. This method requires no random number generation. 
D. COINSS 


COINSS ts designed as a component to be used in conjunction with a map-based 
planning system. The map-based planning system must have certain capabilities and 
provide distinct information to COINSS. First, the map-based planning system must be 
able to provide positions on the map by latitude and longitude. This is required when 
initially creating the search region, as well as when initial areas and search areas are 
created. Second, it must also be able to provide the number of pixels, both horizontally 
and vertically, between two points on the map. This is required when initially creating 
the search region. Finally, it must be able to receive numerical inputs that represent 


probabilities. This is required for initial area probabilities and search area probability. 


COINSS is written in the Java programming language and is platform 
independent. This means that compiled code can be executed on a variety of computer 
systems that have the Java Virtual Machine (JVM), for example, Windows 95, 98, NT, 
and various UNIX platforms. COINSS can be loaded to any computer with the JVM and 
executed immediately. The map-based planning system used to demonstrate COINSS is 
also written in Java. Therefore, the total system can be distributed to execute on any 


computer with a JVM. 


COINSS is written as a Java application, which means that it is executed from a 
command line prompt. COINSS and the map-based system could be converted to a Java 
applet, which would then be executed within a web browser that has aJ VM, for example, 
Netscape Navigator or Internet Explorer. This would be the most convenient method of 
distributing a dynamic view of the search. This would allow searchers, family members, 


and others to have the same view of the search that the Commander has. 


COINSS uses the stochastic representation for target movement, which is 
explained in complete detail in Chapter III. COINSS and a map-based planning system 
are available by contacting Professor Gordon Bradley at NPS, Monterey, CA 93943 or 


email to bradley @nps.navy.mil. 
E. THESIS STRUCTURE 


Chapter If explains the search and rescue process and discusses how the On Scene 


Commander would use a search movement model in a search and rescue operation. 


Chapter III explains the logical search and movement elements of the Coordinated 


Inland area SAR System (COINSS). 
Chapter IV explains the graphical implementation of COINSS. 
Chapter V provides the conclusions and recommendations. 


Appendix A provides the Markov transition matrices for the five classes of targets 


defined in COINSS. 


Appendix B provides a walkthrough of COINSS operating with a map-based 


planning system. 


Appendix C provides the code for SARgrid, and SARhex. This is the complete 


code for COINSS. 





Il. THE SEARCH AND RESCUE PROCESS 
A. INITIAL ACTION 
1. Obtaining an Initial Position 


The initial area is the last known geographical location of the target. The initial 
area can be multiple disjoint or overlapping locations, in which case a probability is 
assigned to each portion sequentially, corresponding to the probability that the target is 
located in that portion. One of three general situations for the initial area information 


usually exists [18]. 


The first is a known position. This position can come from a direct witness, or 


may even come from the target individual. The single position is used as the initial area. 


The second is a known track. In this instance the target’s route of travel is known. 
What is unknown is the position along that route of travel. The complete length of the 


known track is used as the initial area. 


The third is a known general area. This is the least favorable condition, but also 
the most likely. This instance occurs when general information regarding the habits of 
the target is known, such as where he typically hunts, or the part of the lake where she 
typically fishes. Natural boundaries, such as cliffs or rivers, can help to constrain this 


type of initial area. 


In each case, the initial areas are assigned a probability, which is the probability 
that the target is located within that particular geographic area. The combined set of 


locations creates the initial area for the search. 
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2. Defining the Search Region 


The total search region is the geographic area determined to be the most likely to 
contain the target. The bounds of this region need to account for the error in the initial 
position, the movement capabilities of the target, the amount of time the target has had to 


move, and the physical geographic constraints of the area [18]. 


The search region must be large enough to ensure that the target is located within 
it. Conversely, the total search region area should be small to minimize the probability of 
searching areas where the target is not located. Ideally, the total search region should be 
dynamically adjusted, to account for the both the target's movement, and the search 


efforts conducted. 


If available, information about the elevation and terrain within the search region 


can be used to help predict the movement of the target. 
B. SEARCH ACTION 


The search action 1s a looping process. The process evaluates the current 
probability information, identifies the search areas, assigns units to conduct the search, 
executes the search, and then updates the probability information based on the new search 
information. This process is continually repeated, often for a variety of search units, until 
either the target is located or the search is suspended. Because this thesis does not focus 
on the tactical decisions required to create search plans, this section will only discuss the 


role that the outcome of those decisions has on the search movement model. 
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1. Identifying the Search Area 


A search area is a sub-section of the search region that is searched by a particular 
search unit. Search areas are chosen to make best use of the abilities of the search unit in 
order to maximize the probability of detecting the target. Generally, the areas with a 


higher probability of containing the target are searched before the lower probability areas. 
2. Search Execution and Results 


After being assigned a search area, the unit searches that area. Once a search unit 
has completed its search, the unit reports its findings. If the target was found, the search 
can be concluded. If the target was not found, that information is used to update the 
probability map by Bayes Theorem. The probability that the target is in the searched area 
is reduced, and the probability that the target is in any non-searched area is increased. 
The change in probability on the probability map is directly related to the search unit's 
probability of detection of the target, given that the target is located within the area 


searched. 


Upon completion of the search result update, a decision is made whether to 
continue the search or suspend the search. This decision is an easy one when the target 
has been located. When the target has not been located, the decision takes into account 
many factors. These factors include the length of time the search has already expended, 
the probability the target will be located given that it has not been found up to this time, 


and the value of the target. If it is decided to continue, the search cycle is repeated. 
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II. LOGICAL MODEL 
A. PROBABILITY MAP COMPONENTS 


The first layer in the logical model is the probability map. This map presents the 
probability that the target is located at a certain geographic position. This probability 
follows the laws of any standard probability distribution, and computations (such as 
Bayesian updates and using Markov chains) can be applied to it. It is this property which 
allows us to update the map for both time based movement and search result information 


ey. 
1. Elements of the Grid Square 


The logical map can be imagined as a large grid of squares. Each square is an 
individual piece that contains information relating to the complete map. Though the 
information values assigned to each square will vary from square to square, the type of 
information provided will be the same. As such, every square can be represented as an 
individual object, each having its own set of variables. The essential variables for each 
Square are its position, elevation, and a velocity array. The first two of these variables are 
static. These values are set when the model is initialized, and will not change during the 


course of the search. 


The velocity matrix, V, is a four-dimensional array, with the first two dimensions 
representing the location on the probability map, and the second two dimensions being 


speed in the east-west direction and speed in the north-south direction as it would appear 


on acommon geographic map. The units for speed are feet per second. The values of the 
indices indicate the velocity in that direction. For example, Vx 43 represents the 
probability that the target is located in square S, y moving with a velocity of 4 feet per 
second toward the east and 3 feet per second toward the north. Likewise, Vx y_-4.-3 
represents the probability that the target is located in square S, , moving with a velocity 


of 4 feet per second toward the west and 3 feet per second toward the south. 
B. INITIALIZING THE PROBABILITY MAP 


There are two Steps in initializing the logical probability map. The first step is to 
define the entire potential search region. The second step is to identify the potential 


initial areas for the target, as well as the respective probabilities for each area. 


Defining the potential search region is a fairly simple process. Because the 
program uses a grid array, the shape of the region will always be a rectangle. In order to 
define the bounds of this rectangle, the user provides the latitude and longitude of the 
upper left corner of the region and the lower right corner of the region. The map-based 
planning system provides these coordinates to COINSS, along with the length and width 
of the region in terms of graphical pixels. COINSS then determines the search grid array 
size. Currently in COINSS, each grid square is 16 pixels by 16 pixels on the graphical 
display. The search grid array size is then computed by dividing the search region into 


the required number of 16x16 pixel squares to cover the region. 


Defining the initial area for the target is a little more involved. A vector of 


geographic positions, generally denoted by longitude and latitude, is required. These 
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positions will define the vertices for a polygon. This polygon denotes the initial area. 


This polygon does not need to be a regular polygon, and can be non-convex. 


Each polygon is assigned a numerical percentage value, which is the probability 
that the target started in that area. The model takes this probability and distributes it 
equally among the grid squares located within the area. A grid square is considered to be 
located within the area if the center of the square is located within the area, otherwise it is 


considered to be located outside the area. 


Multiple initial areas may exist, and these areas may overlap. Each initial area is 
created in the same manner. It is conceivable that the sum of the probabilities for all the 
initial areas would be less than one. This would occur if there exists a probability that the 


target could have initially started outside of the search region. 
C. UPDATING THE PROBABILITY MAP 


Two events require the update of the probability map. The first is the passage of 
time, which requires an update due to the movement of the target. The second is a report 
from a search team, which requires an update based upon the detection/non-detection of 


the target. 
1. Time-Based Update 


Time based updates adjust the probability distribution on the map due to the 
movement of the target. These updates can occur over any time interval. Smaller time 


intervals, such as one second, come closer to simulating a continuous time movement 
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model. However, larger time steps reduce the computations, and over a large region with 


several hundred movement vectors, the loss of resolution 1s not significant. 


COINSS uses a constant time step, currently set to fifteen minutes. COINSS only 
computes the time-based update when prompted by the user. This allows the user to 


control the pace of the model. The time-based update consists of three major steps: 


1. Moving the velocity array probability elements from their current square to a 


new square, based upon their velocity and the time step increment. 
2. Adjusting for changes in velocity after each probability element movement. 
3. Correcting the velocity array probability elements for elevation. 


The first step of the time-based update is to move the velocity array probability 
elements. This is done by sequentially accessing the velocity array of each square, Vx yi}. 
Each element of the velocity array is moved the direction and distance that is dictated by 


the velocity vector and the corresponding time step increment, T), to a new location, 


S xnew,ynew- 


xnew = x+i*T, 


ynew = y+ j*T, 


The probability at this new location is evenly distributed over a square region, 


centered on the new geographical position (Figure 3-1). The probability is then assigned 
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to each grid square in proportion to the area of the probability square that is located 


within the grid square [6]. 





Figure 3-1. Movement and distribution of velocity array 
probability element 


The next step is to adjust the current velocity of the probability element. The 
velocity of the target is modeled to change according to a Markov chain. To model the 
velocity changes as a Markov chain, each velocity element is recognized as a state. The 
probability of changing velocity from state 7 to state inew, which can be read as changing 
from velocity 1 to velocity inew, is represented by a Markov transition matrix, MMi inew 
[13]. The transition probabilities of this matrix are based on the target type of the search. 
Appendix A discusses the five target types defined in COINSS and provides the Markov 
transition matrix for each type. Figure 3-2 is the transition matrix for target type 1, which 


is a child age 1-6. 


Transition matrix (feet per second velocity) 

4 3 -2 -1 O 1 2 3. 4 
-4 |0.01| 0.06/0.43|0.43]0.06/0.01| 0] 0} _—~O 
-3 | _0/ 0.01/0.11/0.76/0.11]0.01| 0} 0] 0 
-2 |__ 0] 0.01|0.07| 0.6/0.27/0.04/0.01| O|_—O 
-1|__0| 0.01/0.04/0.27| 0.60.07/0.01| _0|_—0. 
O {| _0/0.01/0.11]0.76/0.11/0.01| oo] 
1 | 0] _0/0.01/0.07| 0.6/0.27/0.04/ 0.01| 0 
2 | 0] _0/0.01/0.04/0.27| 0.6/0.07] 0.01| _O} 
3 | O| 0] _0/0.01/0.11/0.76/0.11/ 0.01] 0: 
4 | 0} 0} _0/0.01}0.06/0.43/0.43] 0.06]0.01) 





Figure 3-2. Transition Matrix for a child age 1-6. 
Each row of this matrix sums to 1. 


This matrix is applied individually to both the x and y components of the velocity. 
For example, from Figure 3-2, the probability of transitioning from velocity -2 (ft/sec) to 
-1 (ft/sec) is .6, and is written as MM ..2.;. The probability element, labeled V, i,j, can 
then be distributed to the velocity array of its new square, Vxnew,ynew,inew,jnew, according to 


the following: 


= *k *" 
vee ynew,inew, jnew » Ve yey MM i,inew MM j,jnew 
x'=xnew-T,*i, 


y =ynew-T,*j, 
i,j 


For example, suppose V,,y,-2.2 =.25. Then suppose that by following the example 


of Figure 1, it is computed that 40% of the moved square lies within Sxnew,ynew, OF Vx’y’.i,j 
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= .1, where x’=xnew-T,*1 and y’=ynew-T*j. By using the transition matrix for target 


type | from Figure 3-2, we obtain the resulting distribution shown in Figure 3-3. 


jnew 
-4 -3 -2 -1 0 1 2 3 #4 


inew 





Figure 3-3. Resulting distribution Of Vxnew,ynew,inew,jnew When 40% of 
Vx,y,-2,-2 =-25 1S moved to square Synew,ynew, for target type 1, 
small child age 1-6. The sum of all elements in this matrix is .1. 


Once the probabilities have been moved and updated for all the squares in the 
search region, the final step of the time-based update adjusts the velocity array 
probability elements to account for elevation. For COINSS, elevation can be loaded from 
a CD-ROM that contains Digital Terrain Elevation Data Level 1 (DTED1!) from the 
National Imagery and Mapping Agency. This data provides the elevation in meters and 
is available to the U.S. Government only. The coverage for this data is one measurement 
every three arc seconds, or approximately every 300 feet. If elevation data is not 


available, then the elevation is considered to be O meters for all locations. 


The adjustments for elevation are done by again cycling through all the squares of 


the search region. At each square, the total probability of that square is computed. 
Pinit, , = Ee 
tJ 


oN 


The next step is to cycle through the velocity array V of that square. At each 
probability element, V, ij, of the array, a slope is computed. This slope is computed in 
the direction that the element would move during the next time step (i.e., the slope from 
Sx,y tO Sxnew,ynew , Where xnew = x +1*T, and ynew = y +j*T)). If the slope is too steep for 
the particular target type to easily navigate in that direction, then the probability element, 
Vxy,ij. 1s reduced. If the probability is reduced, it is either reduced to 50% of its previous 
value or to zero, again depending on the severity of the slope and the target type (Figure 


3-4.). These values were determined from judgements of human tendencies. 
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Figure 3-4. Slopes at which the probability of moving is 
reduced either 50% or 100% 
After all elements of the array V in Square S, , have been updated, the total 
probability of the square must be readjusted back to its original value, Pinit, , by 
computing the new total square probability, and re-scaling each element of the velocity 


array. Let V', yi; be the velocity array elements after adjustment for elevation, and let 
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2. Search-Based Update 


Each search conducted in the search region is applied to the model independent of 
any other search. The variables required for a search are the search unit name, the search 
area, and the probability of detection. The search unit name is a label to help identify this 
particular search instance. The search area is a vector of geographic positions, generally 
denoted by longitude and latitude. These positions will define the vertices for a polygon. 
Unlike the search region, the search area is not required to be a rectangle and it may be a 
non-convex polygon. The probability of detection 1s a value between 0 and | that is 


provided by the user. 


Search based updates adjust the probability distribution for all the squares in the 
logical map due to search unit results. These updates use Bayes theorem to account for 
unsuccessful search results. This result is computed by using conditional probability rules 


and Bayes theorem. The following two events are defined: 


D, y: The target is detected at location S, by the search effort allocated to S, , 
Ty: The target's current location is S,, 


Then the probability of detection for the search unit, conditioned on the event that 


the target is located 1n the searched area, is defined as: 


P[Dx.y | Tx,y] 


and the probability of non-detection, conditioned on the event that the target is located in 


the searched area is defined as: 
BIDDERS |= 1b Dei 
where the superscript denotes negation. Likewise: 


P[D,, 





ee a Q, V(x, y) oe (xe y’) 
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Sx,y is determined to be located within the search area if and only if the center of 


square S, y is located within the search area polygon. 


When a Search party reports that they have searched an area and did not find the 
target, the model is updated. Each square in the search region must first be updated to 
show the probability that the target is located in that region and was not detected by the 


search party. This is defined as: 


PIT, y Dy, 1= PLD}, 





eer eee 


Using Bayes theorem, we arrive at a posterior probability that is the probability of 
the target being located at location S,.y, given that it was not detected at location S,,y. 


This is defined as: 


el esl alee eae 


Pilg |= —— oy 
eo alle alee 


It is worth noting that all squares in the search region are being updated, not just 


the squares located within the search area. For all the locations S,:y located outside the 





search area, P[D‘, ,|T,.,.]=1, and therefore P[T,., (1D,.,]=PIT,.,.]. These values are 


used in the above equation to create a new probability distribution for the location of the 


target. 
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IV. VISUAL REPRESENTATION OF THE MODEL 
A. COMPONENT-BASED DESIGN 


COINSS was designed to be a component-based system. COINSS provides the 
algorithms, methods, and logic for executing a SAR mission plan. It is designed to take 
advantage of well-defined reusable components from existing map-based planning 


systems, such as the ability to read latitude and longitude or to provide overlay displays. 


The "Loosely Coupled Components” project designed by Dr. G.H. Bradley and 
Dr. A.H. Buss, both of the Naval Postgraduate School, provides one such map-based 
planning system [4]. Their project creates planning systems by integrating various 
components as required. Components can be maps, algorithms, data overlays, or any 
other item necessary to complete a task. The nature of the system 1s such that the planner 
only needs to assemble the distinct components needed at that time. If a need for a new 
component arises during a planning operation, the planner can obtain or implement that 
component and add it to the planning system without having to restart the entire 
operation. It also allows designers to concentrate on creating algorithms or methods to 
solve problems, and not continually "re-invent the wheel” by developing a graphical 


interface for their specific problem. 
B. BASIC CONSTRUCTION OF THE VIEW 
1. Grid Dimensions 


To present a COINSS probability map on a visual display, a translation from the 


squares on the logical probability map to the pixels on the monitor is required. A 
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standard method of presenting probability maps on a visual display is by using different 
colors or shades of a single color to represent the probabilities. The visual display square 
size needs to be large enough to allow visual discrimination of color and the squares must 
be small enough to provide useful search information (1.e., a search region with only two 
Squares would not provide much meaningful search information). Finally, the logical 
array needs to be small enough to provide movement and search updates in a reasonable 
amount of computation time. COINSS was developed on a PC with a 166 MHZ CPU. A 


16x16 pixel square is a reasonable choice for this computer. 
2. Probability Display 


COINSS uses 10 shades of the color red to discriminate among the different 
quantitative probability values. The red shades are equally spaced, which is achieved by 
manipulating the RGB color value of the pixels. The RGB value ranges from 0 to 255. 
The shades were created by modifying the values in increments of 25 from 5 to 255. Red 
was chosen as the color to use, and the darkest red was achieved with RGB values 


2551010, while the lightest shade was achieved with RGB values 255|225/225 [8]. 


The variable maxSquareProb contains the value of the largest probability in the 


search region. 


maxSquareProb = max{> Veyagir VY 
ij 


Each square is assigned a color, which is based on the probability of the square 
compared to the maxSquareProb. The maxSquareProb is assigned the purest color value 


(2551010). Squares containing a probability that is 90%-99% of maxSquareProb receive 
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the next highest purest value (255/25!25). This continues at 10% intervals down to the 
squares that have probabilities that are between 10%-19% of the maxSquareProb. These 
squares receive the least pure color value (25512251225). The squares with probability 
less than 10% of the maxSquareProb are assigned to be completely transparent, and 


appear as blank squares on the probability map. 


COINSS also makes use of opaqueness. The opaqueness of a cell determines how 
well you can see through the cell on the top layer and view images on the layer(s) below 
it. An opaque value of 255 is completely opaque and does not allow the viewer to see 
any part of the layer below it, while an opaque value of 0 1s completely transparent [8]. 
COINSS currently assigns each colored square in the probability map an opaqueness 
value of 190. This allows the layer immediately below the probability map to be seen, 


with little loss in color discrimination. 
C. INTERACTING WITH THE VIEW 
1. Selecting a Geographic Map 


COINSS is designed to operate in conjunction with a map-based planning system, 
such as the Loosely Coupled Components system created by Bradley and Buss [4]. The 
map-based system is responsible for obtaining a geographic map from any available 
source. The planning system is also responsible for providing an interface to COINSS, 
specifically, a reference to the latitude and longitude of the locations on the map, as well 
as a representation of the scale of the map provided in pixels per unit longitude. The map 


can have any geographic or topological features or symbols. This map is merely a visual 
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reference for the planner, and will have no impact to the logical computations that take 


place in the model. 
2. Defining the Search Region 


Identifying the bounds for the search region 1s the next step. This is a user- 
defined region that is provided by the map-based planning system. The search region is 
required to be a rectangular region, and is defined by two points, the upper left corner and 
the lower right corner. A simple method is to allow a mouse click at the location for the 
upper left corner, followed by another click at the location for the lower right corner. 


These two values are then provided to COINSS to define the search region. 
3. Defining the Initial Areas 


Initial areas are also user-defined areas that are provided by the map-based 
planning system. Again the mouse is used to outline the initial areas. As discussed 
earlier, a vector of geographic positions is required. These positions define the vertices 
for a polygon. The mouse clicks define each vertex, in order. Once the polygon is 
completed, it is assigned a numerical] value, which is the probability that the target started 
in that area. COINSS takes this probability and distributes it equally among the grid 


squares located within the region, and across all the possible velocity vectors. 


This is an iterative process, conducted on each initial area in a sequential manner. 
Each subsequent initial area probability is added to the existing probability map. If two 
or more initial areas overlap, then the probability associated with the overlapped area is 


equal to the sum of the individual probabilities for that square. 
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4. Defining the Search Area 


Search areas are the third case of user-defined areas. These are also provided by 
the map-based system, and are created in the same manner as the initial areas. The 
numerical value assigned to each search area is the probability that the unit searching that 
area will locate the target if the target is in that area. This probability is used to apply 


Baysian updates to the probability map. 


eo 








V. CONCLUSIONS AND RECOMMENDATIONS 


COINSS provides an inland search and rescue model which implements target 
movement during the search execution. Target movement is modeled as a discrete time 
Markov chain, based upon the type of the target. Geographic elevation is also accounted 
for in the movement model by limiting movement in areas with large changes in 


elevation. 


COINSS is a component-based system, designed for implementation with a map- 
based planning system. This design demonstrates that a component-based design can 
reuse components. Allowing different component systems to use the same map-based 
system interface increases the efficiency of the hardware, reduces user training, and 


reduces software cost and redundancy. 


This thesis provides several areas for future research. The movement model is 
very rudimentary. The COINSS Markov movement model is based upon historical data 
that only provides the total distance traveled from initial datum [9][17]. More detailed 


study of target movement rates and tendencies would improve the movement model. 


Integrating additional information regarding the terrain would also improve the 
movement model. Identification of impediments such as heavy brush, rivers, and lakes, 
or easier traveling terrain such as trails, paths, and roads would improve the realistic 


modeling of movement, and create a more realistic probability model. 
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APPENDIX A. MARKOV MOVEMENT MATRICES 


At each movement phase update, the velocity vector of the target changes as 
discussed in Chapter III. To implement the movement model using a Markov chain, a 


discretized version of a Kalman filter was created [22]. 


Using a discrete version of the Ornstein-Uhlenbeck (OV) process [22], velocity is 


modeled as 
Xt = CX, oe Ww, 


where X, is the velocity of the target at time t, W; is an independent identically 
distributed normal random variable with mean O and variance q, and c is aconstant. The 


values for c and q determine the variance of the velocity for the target. Furthermore, 


where A is the time interval at which updates are made, and 1/t is the expected rate (in 
time units) that the target changes velocity [22]. A is an artificial variable which is 
determined by COINSS, while t is a characteristic of the target that is estimated from the 


data. 
The variance of the total movement of the target can be expressed as [19] 


} é 


Var(X (0)) = = (2s*r)[26 —7[3 PAG a) 
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where s” is a target characteristic which identifies the variance of velocity, and 6 is the 
time interval over which the target has moved. By computing the variance for different 
target types from Syrotuck [17] and Hill [9] (Figure A-1.) and estimating values for 6 and 


t, S- can be solved for, and then q can be computed from [22]: 


q=s°(1-c’) 


Target Type Hill Var 
] 


2.15 
30.05 10.62 
22.47 





Figure A-1. Variance of total movement (in miles) for found targets 
in Hill[9] and Syrotuck[17]. 

The data provided by Hill and Syrotuck in Figure A-1 have two shortcomings. 
First, the data is only for found targets. As such, we must assume that all targets move in 
a fashion similar to that of the found targets. Second, the data has no time associated 
with it. The longer a target remains undetected, the larger the variance of his movement. 
However, this is the best data available on lost targets. Therefore, the data was used, 
along with a generic view of human movement patterns to create one model of 


movement. 
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After c and q have been estimated, the Markov transition matrix can be 
constructed using the assumption that velocity changes will vary according to a standard 
normal distribution. Figure A-2 provides a chart of all the values for t, which were 
estimated based upon human tendencies, and a which was computed with an estimated 6 
of 24 hours. Figures A-3 through A-7 are the transition matrices for each target type © 
when A = 15 minutes. The indices for the transition matrices are in feet per second. The 


target types are [9]: 
1: children under age 6, 
2: children age 7-12, 
3: despondents/walkaways, 
4: hunters, 
a: hikers/adults. 


Target Type (sec) S* (ft*/sec*) 
i ia 
1800 


Figure A-2. Computed values for target movement 
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4-3 2 -1_0 1 2 3 4 
0.020.14/0.34/0.34]0.14]0.02| 0/0} 
| 010.04]/0.24/0.44]0.24]0.04| 0} 0] 0 
| __010.03/0.19/0.39]0.29]0.09)0.01| 0] 
|__010.01/0.11/0.31/0.37/0.17/0.03| 0] O} 
| 0} 0/0.04/0.24/0.44/0.2410.04| 0] __O: 
| 0}  0/0.03/0.17/0.37/0.3110.11]0.01| 0} 
| O| _ 0/0.01/0.09/0.29/0.39/0.19]0.03|__0 
| 0} 0} 0/0.040.24)0.4410.2410.04| 0! 
| O} 0} _0/0.02)0.14)0.34]0.34]0.14/0.02: 















Figure A-3. Transition matrix for target type 1 





Figure A-4. Transition matrix for target type 2 
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Figure A-5. Transition matrix for target type 3 


4-3-2 et 2 
S2[oBeloial—o[ of of of o[ 0 
0.04/0.33|0.53| 0.1| Oj O| O| O| O 
-0|0.07|0.43)0.43/0.07| |__| _o| oO 
| Ol 0/0.14/0.48/0.28] 0.1] 0} 0} 0 
| 0] _0/0.02/0.23/0.51/0.23/0.02| Ol. 
| Ol} 0] 0} (0.1/0.28/0.48/0.14; Oo] 0 
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Figure A-6. Transition matrix for target type 4 
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Figure A-7. Transition matrix for target type 5 
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APPENDIX B. COINSS SCENARIO 


This is a walk through of COINSS implemented on a map-based planning 
system. The system used to implement COINSS in this scenario is Special 


Operations Forces/Loosely Coupled Components (SOFLCC) [2]. 





ber a 


Figure B-1. The standard map-based planning system. The latitude and longitude of the 
mouse position is displayed in the lower right of the display. 
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Figure B-2. To add a new overlay to the map, click the “Add a new Overlay” Icon. This 


then provides the option to add a (new) blank layer, or add a (pre-existing) layer from a 
file. 
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Figure B-3. The search region is created by clicking once on the map at the north-west 
corner of the search region, then double-clicking at the south-east corner of the region. 
The map has been faded to provide a clearer view of the search region. 
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Figure B-4. Select the Show Overlay icon to select the search region for initialization. 
Select the Run Algorithm option to execute a SAR algorithm. 
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Figure B-5. Select the Add Algorithm option to dynamically load the SAR algorithms. 
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Figure B-6. Select the SAR .initialRegion to create the search region. 
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Figure B-7. The initialRegion algorithm requires a name as an input parameter. 
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Figure B-8. The system queries for the type of search target, which is contained in a file. 
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Figure B-9. If elevation data is available, the user can import that data into the system. 
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Figure B-10. Add a new overlay to create the initial area. The area is created by single 
clicks of the mouse at each location of the map which is a vertex for the area. Double 
click at the last vertex to complete the area. 
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Figure B-11. Use the Show Overlay icon to select the initial area graph, and run the 
initialArea algorithm. This algorithm requires a name and a probability between 0 and 1 
as parameters. 
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Figure B-12. The probability map updates to show the current probability distribution. 
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Figure B-13. Use the Show Overlay icon to select the search region, and run the 
movement algorithm. The algorithm requires a name and a length of time in minutes. 
The probability map is updated based on the target movement that would occur during 
that time interval. A second movement option ts available, the sleep update. This option 
also requires a name and a length of time in minutes. The time on the probability map is 
updated, but no target movement occurs, because the target is assumed to be sleeping. 
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Figure B-14. Add a new overlay to create the search area. The area is created by single 
clicks of the mouse at each location of the map which is a vertex for the area. Double 
click at the last vertex to complete the area. 
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Figure B-15. Use the Show Overlay icon to select the search area graph, and run the 
searchArea algorithm. This algorithm requires a name and a probability between 0 and 1 
as parameters. The probability map is updated based on the search results. 
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Figure B-16. Double clicking on the probability map displays the current date and time of 
the search. 
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APPENDIX C. COINSS CODE 


This appendix contains the code for SARhex, and SARgrid. In addition, it 
discusses the driver program necessary to interface COINSS with a map-based planning 
system. SARhex is the basic building block for COINSS. SARgrid contains all the 


methods which operate the model. 


The COINSS software is connected to a map-based planning system through a 
driver program that invokes the methods in COINSS. The driver invokes these methods 
by obtaining the necessary information from the map-based planning system. The driver 
then converts the information into the correct parameter form to call the methods in 
COINSS. The driver then invokes the methods using the proper parameters, and provides 


the results back to the map-based system. 


To initialize the search region, the driver instantiates an object “SearchRegion” of 
class SARgrid. The constructor requires the latitude and longitude of the northwest 
comer and the southeast corner, each being double variables with units degree. The 
constructor also requires the height and width of this region on the screen in pixels, each 


being int variables. 


To specify the initial area(s) the method probInitialDatum is invoked on the object 
SearchRegion. The first parameter is a two-dimensional array of type int. Each column 
of the two-row array is a vertex of the polygon that defines one initial area. Each column 
is the location, in grid squares (x,y), of a point. The second parameter is a double 


variable that is the probability (0.0<j < 1.0) that the target started in this area. 


2) 


To specify the search area(s), the method probInfoUpdate is invoked on the object 
SearchRegion. The parameters are the same as the method probInitialDatum, where the 
double variable is the probability the search team will find the target, given he is located 


in the given search area. 


To specify movement, the method probTimeUpdate is invoked on the object 
SearchRegion. This method requires no parameters. It will update the probability map 
for movement of one time step, based upon the constant time step defined by SARgrid 


(currently 15 minutes). 


To specify time passing during the night, the method updateTime is invoked on 
the object SearchRegion. The only parameter is a double variable. This double is the 


number of seconds that have elapsed while the target sleeps. 


To present the probability map, the method getProbMapIcon is invoked on the 
object SearchRegion. This method requires no parameters. It returns the current 


probability map as an Icon variable. 
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